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Abstract 

Motivated by modern applications like image processing and wireless sensor networks, we 
consider a variation of the Kepler Conjecture. Given any infinite set of unit balls covering the 
whole space, we want to know the optimal (lim sup) density of the volume which is covered by 
exactly one ball (i.e., the maximum such density over all unit ball covers, called the optimal 
1-density and denoted as 6d, where d is the dimension of the Euclidean space). We prove that 
in 2D the optimal 1-density S2 = {3^/3 — tt)/tt ^ 0.6539, which is achieved through a regular 
hexagonal covering. In 3D, the problem is widely open and we present a Dodecahedral Cover 
Conjecture which states that the optimal 1-density in 3D, 63, is bounded from above by the 
1-density of a unit ball whose Voronoi polyhedron is a regular dodecahedron of circum-radius 
one (determined by twelve extra unit balls). We show numerically that this 1-density S^^dc) 
« 0.315. 
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1 Introduction 



The Kepler conjecture asserts that the density of a packing of unit balls in three dimensions (3D) 
is no more than vr/VTs ~ 0.74048. This is known as the oldest (unsolved?) problem in discrete 
geometry and is listed as part of the Hilbert's 18th problem. About ten years ago, Hales published 
a series of papers concerning the proof of the Kepler conjecture, using a vast amount of computer 
verification [7]. While many believe that the conjecture has been proven, some others believe it 
has not (as the referees for Annals of Mathematics claimed that "The news from the referees is 
bad... They have not been able to certify the correctness of the proof, and will not be able to 
certify it in the future" ^12]). In fact, some of Hales' computer programs (total 3GB!) need to re- 
checked, which will take another 20 years ( jhttp : //en. wikipedia . org/ wiki/Kepler_con j ecture ) . 

In this paper, we make no attempt in re-proving the Kepler conjecture. In fact, this research is 
initiated by looking at the Kepler conjecture from another perspective: What is the origin of the 
Kepler conjecture? Application! It is known that T. Hariot first investigated the Kepler conjecture 
in 1591 while working for Sir Walter Raleigh, who asked him to determine formulas for the num- 
ber of cannonballs regularly packed. Later on Hariot communicated the problem to Kepler. This 
problem on unit ball packing is known as the Kepler conjecture since [7]. 

Looking at the story, it is not difficult for one to realize that applications do drive theoretical 
research. In the cannonball packing problem, obviously cannonballs cannot overlap. But in some 
modern applications, the situation could be different. Example 1: in wireless sensor networks, 
many homogeneous sensors, whose sensing regions are modeled as unit disks in 2D, are placed 
in a region R (we can assume that the diameter and width of R are much bigger than one). As 
sensors have limited battery time, to extend network coverage lifetime (which is the time that the 
region is completely covered) one needs to schedule the sensors, i.e., turn some on and some off. 
The requirement is that all the time when the network is alive one must make sure that every 
point in R is covered by some sensor. The abstract version of this problem has been shown to be 
NP-complete [3], but the complexity of the geometric version of the problem is still open (though 
it seems possible to extend the technique in [3J to obtain a PTAS). In any case, in this wireless 
sensor network problem clearly the sensing regions can overlap. 

Another related example is regarding the conversion of continuous-tone images into binary 
images for printers, formally known as the digital-halftoning technique. Given a continuous-tone 
image, one can compute the spatial frequency distribution by applying the Laplacian or Sobel 
differential operators. Then we have a grid of points of possibly slightly different radii (or disks 
centered at grid points), and the problem is to compute a subset of disks which maximize the area 
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of the regions covered by exactly one disk. The complexity of the problem is still open, a factor-5.83 
approximation was given in [1] . When the ratio of the radius of the largest disk over the radius of 
the smallest disk is a constant, then the problem admits a PTAS [3]. 




Figure 1. A regular hexagonal covering. 

We feel that the reason why we still could not completely solve these kinds of geometric problems 
is that some fundamental question regarding disk covering, very surprisingly, is still untouched. We 
cover one such problem in this paper; namely, the 1-density of any unit disk cover which covers 
all the points in the plane (2D) Qj. We prove a theorem which is corresponding to Thue's theorem 
on the density of unit disk packing; loosely speaking, the regular hexagonal covering achieves the 
optimal 1-density at 0.6539 (Figure 1). (A regular hexagonal covering is the configuration where 
the Voronoi polygon of any center of the unit disks is a congruent regular hexagon. Readers are 
referred to [TT] for concepts and properties regarding Voronoi diagrams.) 

We then consider the problem in 3D space, i.e., given any unit ball cover covering all the points 
in space, decide the corresponding optimal 1-density. To some extent, this is in fact similar to the 
Kepler conjecture, which is on unit ball packing. But the status of this 3D problem is widely open. 
Our numerical calculation with computer programs indicates that the 1-density of a unit ball whose 
Voronoi polyhedron is a regular dodecahedron with circum-radius one (determined by twelve extra 
unit balls) is around 0.315. We conjecture that this is the upper bound of the optimal 1-density 
^3, which is also in parallel with the Dodecahedral Conjecture by L. Fejes Toth |13l [8] regarding 
the Kepler conjecture. We will call this the Dodecahedral Cover Conjecture henceforth. 

This paper is organized as follows. In Section 2, we present necessary definitions and a proof 

^Note that the density for finite unit disk (ball) covering has long been studied [2- But this traditional density 
is not related to the 1-density we will discuss in this paper. The traditional density for unit ball covering refers to 
the maximum ratio of the volume of some convex region covered by the balls over the volume of the whole region 
covered by the balls. 
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for the 2D case. In Section 3, we show how the 1-density in the Dodecahedral Cover Conjecture 
can be computed (estimated) empirically. In Section 4, we discuss our conjecture and leave some 
other open problems. 



2 Result in 2D 

We first make some definitions for any d-dimension. 

Definition 2.1 A unit ball cover in d- dimensional space is an infinite set of unit balls which cover 
every point in the d-space. 

Definition 2.2 Given any unit hall cover C , let the subset of the halls in C which covers a larger 
hall B with a fixed center and a radius r > 1 be C{B,r). Let Ai{C{B ^r)) be the volume of the 
parts in B which is covered by exactly one hall in C{B,r). Let A{B) he the volume of B. Then, 
the 1-density ofC, 6(i{C), is defined as 

5d{C)= hm sup J— . 

Among all possible unit ball covers, the one achieving the maximum 1-density is denoted as C* and 
the corresponding optimal 1-density is denoted as 5^- 

For the 2D case, a 'ball' is usually called a disk and 'space' is typically called a plane. The main 
theorem for the 2D case is based on the following lemmas. 

Lemma 2.1 Let P = {pi,P2)P3j •••} be an infinite set of points in the 2D plane such that P lie 
on an infinite (or unbounded) number of lines. Let YTi(P) be the Voronoi diagram of P such that 
every point q in the plane is covered by some unit disk centered at pj. Then the Voronoi polygon of 
Pi, VP (pi), is bounded for alii. 

Proof. If VP{pi) is unbounded, then pi is on the convex hull of P = {pi,P2,P3, ■■■}■ Then some 
points in VP{pi) is not covered by any unit disk centered at pj for any j. A contradiction. □ 
Let D = {Di,D2,Ds, ...} be an infinite number of disks of a uniform radius, with each disk Di 
centered at Oj. Let VD{D) be the Voronoi diagram of the center Oj's for the disks in D. For any 
Oj, let Vi be the Voronoi vertex in the Voronoi polygon VP{oi) which is the farthest from Oj. Define 
^max = supj. d(ojfc, Ujfc). We have the following lemma. 

Lemma 2.2 D = {Di\i = 1,2,3, ...} is a unit disk cover if and only if Vi is covered by Di, for all 
i. 



3 



Figure 2. Illustration for the proof of Lemma [2T2l 



Proof. We first show the "if part. Without loss of generality, let r^ax = 1 > '"i = d{oi,Vi), for any 
i. By the definition of Vi, a disk !)■ centered at o, with radius covers VP{oi). As Dj is centered 
at Oi and with a radius r^ax = ^ > ri, Di also covers VP{oi) (Figure 2). So every point in the plane 
is covered by some unit disk in D. D is hence a unit disk cover. 

We now show the "only if part. If D is a unit disk cover, then by definition, every point in 
the plane is covered by some disk Dk- Following Lemma l2.ll all the Voronoi polygons VP{oi) in 
VD{D) are bounded, for all i. Without loss of generality, let Tmax = 1 = d{oi,Vi), for some i. By 
definition, Vi is the farthest Voronoi vertex in VP{oi) from Oj. Therefore the disk Di covers Vi (and 
VP{oi)). □ 

Theorem 2.1 S2 = (3^3 - 7r)/7r 0.6539. 

Proof. Let C = {Ci, C2, C3, ...} be any unit disk cover, with each disk Cj centered at Oj. Without 
loss of generality, let rmax = 1 = d{oi,Vi), where Vi is the farther Voronoi vertex on VP{oi) 
from Oi. By the definition of a Voronoi vertex, Vi is shared by at least three Voronoi polygons 
1^. Assume that one of the other Voronoi polygons containing Vi as a vertex is VP{oj), we have 
d{oi,Vi) = d{oj,Vi) = 1 (notice that Vi is the farthest Voronoi vertex on VP{oj) from oj, or we can 
imagine that Vj = Vi). Following Lemma \2.2\ to cover the whole plane we must cover Vi using a 
unit disk centered at Oj. Let d{oi, oj) = 2x. Denote the sector in Cj bounded by rays Oivl and Ojoj 
as (viOiOj) (Figure 3). Obviously, to maximize the 1-density of sector (viOiOj), it is reasonable to 
assume that no disk Ck intersects the region in sector (viOiOj) which is out of the disk Cj. This 
part of the area in (viOiOj) which is only covered by Ci, is 
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By a simple calculation (taking the derivative over x then solving a simple equation), this function 
is maximized at when we set x = which implies that the angle between Oivl and Ojoj is 



We do not assume that no three centers are collinear, as long as all the centers are on an infinite number of lines. 
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7r/6 (Figure 3). Correspondingly, the optimal 1-density in the sector (viOiOj) is 

3\/3 — vr vr 3\/3 — vr 
^ 12 ^'^W " ^ ' 

which is roughly 0.6539. By symmetry, the 1-density for Cj is optimal when Cj and the five other 
unit disks intersect Ci in a regular hexagonal setting (e.g., as in Figure 1). Clearly this configura- 
tion can be arranged at all unit disks. 




Therefore, 

3^/3 -vr 



0.6539. 



vr 

□ 

We comment that the above optimal regular hexagonal covering is similar to Thue's theorem 
[14:\ I15j . which states that the densest density of unit disk packing is achieved with a regular 
hexagonal packing, which admits a density at ~ 0.9068. 



3 Result in 3D 

For the 3D problem, the most relevant application seems to be placing sensors (or sonars) in the 
ocean to monitor the temperature change or the appearance of submarines, etc. It turns out that 
for the 3D problem, it might be difficult to apply the same technique as in Theorem l2.1l (as the 3D 
counterpart of the regular hexagons, i.e., regular dodecahedra, cannot even tile the whole space). 
Another reason is the difficulty in computing the volume of a non-convex body (i.e., it might not 
even be a polyhedron) analytically. This is very similar to the Kepler conjecture. However, different 
from the Kepler conjecture, the status on computing/estimating the optimal 1-density ^3 is widely 
open. 

We conjecture that S3 is bounded from above by the 1-density of a unit ball whose Voronoi poly- 
hedron is a regular dodecahedron with circum-radius one (determined by twelve extra unit balls). 
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We denote this 1-density as Ss^dc). As we do not know how to calculate S3{dc) analytically (which 
is related to compute/estimate the volume of a non-convex body), we first use a random sampling 
technique to estimate 53{dc). Incidentally, while there have been a lot of theoretical research on 
estimating the volume of a convex body, estimating/approximating the volume of a non-convex 
body has only been practically investigated recently by the author of this paper and his colleagues 
(though without a theoretical guarantee) [lOJ. (The problem of estimating/approximating the vol- 
ume of a convex body, with some theoretical guarantee, has been well researched [9]. However, for 
non-convex geometric body, the problem is still open.) We will sketch some details of the method 
(called Generate-and-Probe) below and apply it to obtain an approximate 63{dc). We comment 
that Generate-and-Probe is a very general method, empirically it works even when no description of 
the body is given — as long as there is an oracle which tells whether a point is inside the body or not. 



Following Figure 4, the unit ball Bi centered at Oj is a circumscribing ball for its Voronoi 
polyhedron VT(oj) — a regular dodecahedron. Computed by a computer program, the edge length 
a of VP{oi) is a=0. 763934. The distance from Oj to one of the faces of VP{oi) is H=0. 760071. 
Define R as the radius of the circumscribing circle of a face of VP{oi), which is a pentagon. Then 
R=0. 649841. So the ratio ai of the volume of FP(oj) to the volume of Bi is 

12 * (0.5 * a * R * sin(3. 14159 * (54.0/180.0)) * 5) * H * /3 



0.728762. 



4* 3.14159/3 

To compute the 1-density of 53{dc), all we need to do now is to compute the ratio 02 of the volume 
inside a dodecahedron which is covered by only one ball to the volume of the dodecahedron. Then 
S3{dc)=aia2- 




(II) 



(III) 



Figure 4. Illustration for the computation of SQ{dc). 



By symmetry, we consider two unit balls centered at 0= (0 , , 0) and pQ= (0,0, 2H) = (0,0,1. 52014) 
respectively. The intersection of the two balls, at height z = E = 0.760071, contains a regular 
pentagon (face) of the regular dodecahedron corresponding to the unit ball centered at pQ. Let 
j»2= (0,0, 0.760071) be the center of this pentagon. Again, by symmetry, we only need to consider 
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a tetrahedron formed by po,p2 and two other neighboring vertices of the pentagon at height z = 
H. We choose them as pi= (0,0. 649841, 0.760071), p3= (0.618035, 0.200812, 0.760071). Then 
the volume of the tetrahedron T = <>{poPiP2P3) formed by Po,Pi,P2,P3, vol(T), is analyticahy 
computed as 

vol(T) = 0.050877. 

The part of T, which is out of the ball centered at o= (0,0,0), or + + = 1, is what we 
need to compute. We denote this part as S and its volume as vol(S). In Figure 4, we show the 
bottom-up view of the regular dodecahedron centered at po, in (I); the side view along the y-axis 
in (II); and the 3D view in (III). 

The volume of S, which is not even a polyhedron, is hard to compute. We first use the Generate- 
and-Probe method of Liu, Zhang and Zhu [lOj. The method starts by generating n random points 
in S, and then computing a small probing ball (the radius of the probing ball is selected so that 
the probability that the ball contains at least some points in S goes to 1 when n — > +oo). Then we 
use the probing ball to probe as many parts of S as possible and count the average ratio TZ of the 
number of random points enclosed by the probing balls over the volume of all the probing balls (of 
course, over a decent number of tries). Eventually, the volume of S is simply n/TZ. We summarize 
our empirical results over three cases with n=80,000, n=140,000 and n=200,000 respectively. The 
average volume of vol(S) is obtained over 100 tries of the Generate-and-Probe method and the 
standard deviation is also given. 



n 


80,000 


140,000 


200,000 


avg(vol(S)) 


0.0219066 


0.0219515 


0.0221303 


cj(vol(S)) 


1.66214e-07 


6.56446e-08 


6.14723e-08 



Table 1. The average of vol(S) and its standard deviation o"(vol(S)), obtained by 

Generate-and-Probe, over 100 tries. 

It remains to verify that the above result by Generate-and-Probe is accurate, as there is no 
theoretical guarantee for Generate-and-Probe. We try to verify this, again empirically. Luckily, 
the boundary of S can be specified explicitly. Let the tetrahedron which is obtained by cutting T 
with a halfspace z > 1 be T' (the plane z = 1 intersects the segment opo at p4=( 0,0,1)). Then 
analytically, we can compute the volume of T', vol(T'). It turns out that 

vol(T') = 0.0163051 < vol(S). 

The volume of the tetrahedron oipoPiPsP^) is 0.0348169, which is an upper bound of vol(S). To 
obtain a closer upper bound for vol(S) , on the surface of x'^ -\- y'^ + z'^ = 1 inside the tetrahedron 
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o{poPiP3P4) (note that pi,P3,P4 are the three corners of this spherical patch), we generate a set M 
of m points. We then compute the 3D convex hull of M U {pi,P3,P4}- We use the upper hull of 
CH{M (J {pi, ps, P4}) to approximate the patch. The approximate volume of vol (S), which we still 
denote as vol(S), is computed by summarizing the volume of all the tetrahedra <>{pouvw) where 
A{uvw) is a face on the upper hull of CH{M U {pi,P3,P4})- The results are summarized in Table 
2. 



m 


avg(vol(S)) 


cj(vol(S) ) 


max(vol(S) j 


2,000 


0.0219467 


6.52514e-ll 


0.0219879 


4,000 


0.0219645 


1.91297e-ll 


0.0219956 


6,000 


0.0219810 


8.92527e-12 


0.0219994 


8,000 


0.0219843 


5.62743e-12 


0.0220002 


10,000 


0.0219873 


4.46376e-12 


0.0220016 


12,000 


0.0219914 


3.10292e-12 


0.0220046 


14,000 


0.0219931 


2.59317e-12 


0.0220046 


16,000 


0.0219947 


2.63297e-12 


0.0220096 


18,000 


0.0219966 


1.48441e-12 


0.0220060 


20,000 


0.0219970 


2.07155e-12 


0.0220060 



Table 2. The average of vol(S), its standard deviation cr(vol(S)) and its maximum 
max(vol(S)), obtained by the 3D Convex Hull method, over 100 tries. 

Prom Table 2, it can also be seen that the maximum of vol(S), over all of our 1000 tries, has 
never been greater than 0.0220096. Therefore, with these rigorous computation and verification, 
it is safe to say that vol(S) 0.022 hence 02 ~ 0.432414 and d3{dc)=aia2 ~ 0.315. The Convex 
Hull method has much smaller deviations, and is probably more accurate in this case, which is un- 
derstandable. As we mentioned earlier, the Generate-and-Probe is a more general (interior-based) 
approach. Even when no exact description of the body is given, it still works (at least empirically). 
We hence have the following claim. 

Claim. Ssidc) « 0.315. 

4 Closing Remarks 

The Kepler conjecture on the 3D unit ball packing has been bothering mathematicians for about 
400 years and a long story has been developed on it. In this paper, based on modern applications in 
image processing and wireless sensor networks, we consider the problem of computing the optimal 1- 
density of a unit ball cover which is a variation of the Kepler conjecture. We prove a theorem for the 
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2D case that this optimal 1-density in 2D is 62 = (3\/3 — 7r)/7r ~ 0.6539. The 3D problem is widely 
open, though we think that the optimal 1-density is bounded from above by the 1-density of a unit 
ball whose Voronoi polyhedron is a regular dodecahedral of circum-radius one (i.e., this Voronoi 
polyhedron is determined by twelve extra unit balls). We denote this 1-density as Ssidc) and we 
show numerically that 6s{dc) ~ 0.315. We conclude this paper with the following formal conjecture. 

The Dodecahedral Cover Conjecture (for unit ball covering): The optimal 1-density in 3D, 
^3, is bounded from above by 63{dc) « 0.315. 

It is also natural to say that when the dimension d is large, the optimal 1-density in d-dimension 
goes to zero, i.e., lim(i_»+oo dd 0. 
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